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We analyze a model of quantum nets and show it has non-abelian topological order of doubled Fibonacci 
type. The ground state has the same topological behavior as that of the corresponding string-net model, but our 
Hamiltonian has less complicated interactions, and the excitations are dynamical, not fixed. This Hamiltonian 
includes terms acting on the spins around a face, around a vertex, and special "Jones-Wenzl" terms that serve to 
couple long loops together. We provide strong evidence for a gap by exact diagonalization, completing the list 
of ingredients necessary for topological order. 
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Introduction- The theoretical and experimental search for 
condensed-matter systems with fractionalized excitations has 
attracted considerable attention in recent years. One reason is 
the promise of protection against errors in a topological quan- 
tum computer 1 1 , 2|. Another is more fundamental: fraction- 
alization provides a dramatic example of how rich the emer- 
gent long-distance behavior of a many-body system can be. 

Finding magnetic systems with fractionalized excitations, 
in particular those with the non-abelian statistics necessary for 
topological quantum computation, is a difficult but important 
problem. Major progress resulted from the introduction of the 
"toric code" or "quantum-double" models, exactly solvable 
models based on any discrete group (I]. Even more types of 
fractionalized excitations and their anyonic braiding occur in 
the "string-net" models |3|. These Hamiltonians are sums of 
commuting projectors, each of which annihilates the ground 
state. The terms must be highly fine-tuned so that they com- 
mute: for example, the string-net Hamiltonian contains inter- 
actions among 12 spins. Since each can be diagonalized indi- 
vidually, the full gapped spectrum of the Hamiltonian is easily 
computed. The excitations are effectively immobile defects 
and so are non-dynamical. Even though the solvable mod- 
els are obviously quite special, the presence of the gap means 
they exemplify a particular phase, not just a special point. In- 
deed, rigorous proofs U |5l show that the topological order 
necessary for fractionalized excitations persists when perturb- 
ing around the exactly solvable points. 

In this paper we focus on a particularly important phase, 
that with "Fibonacci" topological order This is the simplest 
kind of topological order that is not only non-abelian but uni- 
versal; i.e. any unitary operator on the Hilbert space can be 
approximated to arbitrary precision by braiding the quasipar- 
ticles. We find a Hamiltonian with an exact ground state that is 
simpler than that of the exactly solvable string-net Fibonacci 
model 1 3 1, and which further illuminates the origin of this re- 
sulting topological order We expand upon Ref. 6 by including 
terms in the Hamiltonian with the essential property of cou- 
pling loops in different winding sectors. By extensive numer- 
ical exact diagonalization, we then provide strong evidence 
that the resulting Hamiltonian is indeed gapped. Coupled with 
the work of Refs.|6j|7jon the ground state, our work provides 



an explicit Hamiltonian with doubled Fibonacci topological 
order and dynamical excitations. 

Quantum loops and nets- The ground state of a system 
with fractionalized excitations must satisfy special properties. 
In a magnetic model, excitations are typically attached by 
a defect line or "string". The string itself costs no energy, 
because away from the quasiparticles, it is locally indistin- 
guishable from the ground state itself. The existence of the 
string does however result in anyonic statistics, because ro- 
tating one quasiparticle around another causes the strings to 
interact, even when the quasiparticles are far apart. This pic- 
ture is elegantly illustrated in quantum loop models |8 1, where 
the appropriate topological properties are built in directly. The 
Hilbert space is spanned by non-intersecting string configura- 
tions, and the Hamiltonian includes an energy penalty for each 
"string end". Thus the ground state can only be comprised of 
states containing closed loops. Any state with a string end 
then corresponds to a quasiparticle excitation, and the de- 
fect lines are the strings attaching the quasiparticles. There 
are thus for example three distinct configurations with four 
quasiparticles, corresponding to how the strings attach them. 
Since away from a quasiparticle a string with ends is indeed 
indistinguishable from a closed loop, these three configura- 
tions should have the same energy. This degeneracy opens up 
the possibility of non-abelian braiding, because braiding can 
cause a transition between the different configurations. 

Unfortunately, subsequent work |9-11| showed that non- 
abelian anyons do not occur in quantum loop models in their 
simplest formulation. A way of salvaging the idea is to mod- 
ify the inner product to one with more desirable topological 
properties IS). The end result is that when one transforms 
a quantum loop model with the appropriate (non-orthogonal) 
inner product to an orthonormal basis, the excitations are nat- 
urally described in terms of quantum nets, configurations that 
allow branching, or vertices with three or more strings present. 
Thus quite nicely, one ends up at the end with the same sort 
of Hilbert space as in the string-net models ! 

The vertex terms in the Hamiltonian- Our Hamiltonian is 
of "Rokhsar-Kivelson" [12 | or "frustration-free" type. It is the 
sum over (not necessarily commuting) projectors, such that 
there exists at least one state annihilated by each of the projec- 
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tors. Such a state is necessarily a ground state, and all ground 
states must be annihilated by all projectors. Quasiparticle ex- 
citations then necessarily involve configurations not annihi- 
lated by at least one of these projectors; the lowest-energy 
excitations are those dominated by configurations where only 
a few of the projectors acting on it do not annihilate it. 

We study models with a two-state orthonormal quantum 
system |0);, at each link I of some two-dimensional lat- 
tice. The basis elements of the Hilbert space are then the ten- 
sor product of the states for each link. A convenient way of 
picturing each basis state is to draw a line on each link I when 
the state is present, while leaving the link empty for |0);. 
This way, each basis state is pictured by some geometric ob- 
ject; the strings are formed from links with jl);. 

A quantum net Hamiltonian includes an energy penalty for 
"string ends", to ensure none appear in the ground state fl\ [3] 
|8J. This is done by adding a term for each vertex v projecting 
on to any configuration with a string end, i.e. a configuration 
where there is one link / touching v with the state and the 
remainder of links having states with |0)/. Let the projector Pi 
actasPi(- • • |0)/ (g) •••) = (••• (g) \0)i g) • • • ) and • • » 
(g • • • ) = . Thus on the square lattice, with wi, f2, ^3, V4 
the links touching v, the operator 

annihilates any basis state without a string end at v, and gives 
1 acting on any basis state with one. Adding iJ„ to a 
Hamiltonian comprised of sums of projectors therefore for- 
bids string ends in a zero-energy ground state. The zero- 
energy state |\&) annihilated by J^v^v ^'^^ ^^^^ written 
as a sum over "nets" N, geometric objects that have no ends: 
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Other terms in the Hamiltonian then will determine the values 
of the weights w{N) = (A^l^*) in the ground state. 

A ground state with non-abelian topological order- The 
weights w{N) of the configurations in the ground state ^ 
must satisfy certain properties to ensure topological order 
One property is that the number of ground states depends on 
the genus of the surface on which the model is defined. For ex- 
ample, in the toric code HI and in the triangular-lattice quan- 
tum dimer model fT3], there is only one ground state when 
space is topologically a sphere, but there are four when space 
is a torus. A ground state of the form (j2]) is very natural for 
giving this behavior, because geometric objects such as loops 
and nets can have winding numbers around the cycles of a 
surface with non-trivial topology. When the local terms in the 
Hamiltonian do not change the winding numbers, there is a 
ground state for each winding number. 

A gap in the spectrum is also highly desirable. While of 
course a gap is a property of the Hamiltonian, the ground state 
of any local Hamiltonian with a gap in two dimensions satis- 
fies an important constraint: all expectation values of local op- 
erators in such a ground state must decay exponentially as the 



operators are moved apart fT4\. On the flip side, however, to 
have deconfined anyons, the expectation value of some non- 
local operators must decay algebraically; geometric objects 
like loops or nets should have long-range correlations. Ex- 
ponential decay here would confine anyons, because it would 
imply a vanishing probability for two would-be fractionalized 
excitations attached by a string to be far apart. A familiar 
classical model with analogous behavior is percolation: at its 
threshold the probability that two points are on the same clus- 
ter decays algebraically, while by construction all correlators 
of local quantities have zero correlation length. 

To have abelian topological order, a ground state satisfying 
these properties is typically sufficient. To have fractionalized 
excitations with non-abelian statistics, however, much more 
structure need be present. Obviously, there also must be de- 
generacies in the excited-state spectrum so that the exchange 
of quasiparticles can cause a transition to another state. How- 
ever, this is not enough either ll9 UTTIl : the inner product must 
also have the appropriate topological properties |6|. 

In Ref.[l5] it was proposed that non-abelian topological or- 
der occurs when the topological part of ground-state weight 
w{N) is given by a chromatic polynomial. The chromatic 
polynomial XffiQ) easiest to understand by treating the 
strings in the net N as borders separating countries. The dual 
graph N corresponds to a vertex for each country and an edge 
connecting each pair of countries sharing a border For Q in- 
teger, Xn{Q) is the number of ways of coloring each country 
with Q colors such that adjacent countries (i.e. two adjacent 
vertices on N) are colored differently. It is a polynomial in Q 
that can be evaluated for Q non-integer as well, and by def- 
inition is a topological invariant. Any string end results in 
Xj^(Q) = because the dual has a vertex attached to itself. 

We focus on the simplest example of universal non- Abelian 
topological order, the "Fibonacci" case Q = (p"^ = (j> + 1, where 
(f) — {1 + \/5)/2, the golden ratio. The ground-state weight 
w{N) of any net model with (doubled) Fibonacci order, in- 
cluding the string-net model, must necessarily involve the 
chromatic polynomial [161 . This topological order and the 
resulting excitations have been discussed in depth in llT6l[T7l . 
The degrees of freedom in the Fibonacci string-net model lO 
are a two-state system |0);, on the links of the honey- 
comb lattice, and the unnormalized ground state | Vt) in pi is 
summed over all net configurations without ends with 1161 
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where is the number of trivalent vertices in the net N (i.e. 
those with all three neighboring links in the state |1)). To 
find a Hamiltonian simpler than that of the string-net model, 
and which has dynamical excitations, we study a ground state 
slightly different than (pj. Following 16], we take 



w{N) = <j) 



(4) 



where the "length" Ljy of the net is the number of links 
it covers, i.e. the number of states Since the weights 

( 3]4 1 are identical in their topological properties, it is natural 



3 



to expect that the latter results in the same doubled-Fibonacci 
topological order as the string-net ground state (|3]l. Numerical 
and analytical arguments indicate that indeed such a ground 
state has all the desired properties described above IS) IT). 

The face terms in the Hamiltonian- The remainder of this 
paper is devoted to finding a gapped Hamiltonian annihilating 
the ground state \9) in (j2]) with weights (j4]). An advantage 
of using the weights Q is that they obey a "quantum self- 
duality": they take the same form when rewritten on the dual 
lattice 16]. The results in projectors annihilating this ground 
state involving only the links around a face of the lattice, as 
opposed to the 12-spin interaction in the string-net model |[3|. 

Seeing this requires making a change of basis from |0)(, 
on each link to another basis |0)(, via the unitary matrix 
F = defined by 
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Not coincidentally, this matrix is also the fusion matrix for 
four Fibonacci anyons. The quantum self-duality stems from 
interpreting this new basis as describing nets on the dual lat- 
tice. The links of the dual lattice are in one-to-one correspon- 
dence with those of the original lattice, so we are free to de- 
fine such a net D by drawing a string on the dual link when 
the state is present, and leaving the dual link empty when 
|0); is present. The remarkable fact is that when the ground 
state 1^1') is rewritten in the new basis, one obtains li6. .18J 
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D 



w{D) = (Dl*) = a 



(6) 
(7) 



where ijj is the length of the net D, i.e. the number of dual 
links with while a is an unimportant constant. Compar- 
ing (6]7 I with (2]4 1 shows the weighting of the nets of the dual 
lattice is completely equivalent to the weighting on the origi- 
nal lattice, hence the quantum self-duality. 

In particular, (|7]i means that any configuration D with string 
ends has {^{D) = and does not appear in the ground state. 
Thus a projector Hf for vertices on the dual lattice, analogous 
to on the original lattice, will also annihilate this ground 
state. The links /i, /2, . . . touching a vertex on the dual lattice 
correspond to Unks around a face on the original lattice. Thus 
defining the projector Pi = FPiF onto |0)/ means that for 
each original face on the square lattice, the projector 

Hf = a-PMAPu + i^-Ph)PhPhPu (8) 

+ (1 - PfJPhPhPh + (1 - PfJPhPhPfs 

annihilates |^). This operator is non-diagonal in the original 
net basis and so couples different net configurations. 

The Hamiltonian H^f = + J2f Hf thus annihilates 

the desired ground state \'^), and only has interactions around 
each vertex and face, so e.g. only has four-spin interactions 
for the square lattice. The last remaining thing to check is if 
it has a gap. Our numerical results unfortunately indicate that 



there is not a gap for this Hamiltonian, see fig. |5] Another 
(probably related) problem is that Hyf does not have the right 
number (four) of ground states for doubled Fibonacci topolog- 
ical order on the torus. The number grows with the size of the 
torus, which follows by rewriting \'^) in a (non-orthonormal) 
loop basis; the unwanted ground states correspond to loops 
wrapping around the cycles of the torus |6|. 

Coupling wrapping loops via Jones-Wenzl terms — A 
virtue of writing the ground-state weights in terms of topo- 
logical quantities like the chromatic polynomial is that one 
can find many local projectors annihilating it |15|. However, 
many of these will not change the number of ground states 
on a torus. However, "Jones-Wenzl" projectors ||T9| added to 
the Hamiltonian not only couple the long loops and so remove 
the unwanted ground states, but can give a gap JSl. For exam- 
ple, the toric code Hamiltonian can be rewritten in the form 
i?vf plus Jones-Wenzl type terms appropriate for Q = 2 |6|. 
We here derive these terms for the Fibonacci ground state Q, 
and provide strong evidence using numerical exact diagonal- 
ization that the spectrum now includes a gap. 





FIG. 1 : The Jones-Wenzl identity for Q — (fi 

The Jones-Wenzl projectors for knot polynomials re- 
sult in an identity for each chromatic polynomial at 
cos~^{^/Q/2) /tt rational GH. The identity for Q = (t>'^ is l|20ll 



(9) 



with the three nets t,E,I involved illustrated in fig. [T] This 
identity is true locally, meaning it holds for any portion of a 
net. Consider the three nets \t), \E), and |/) on the square lat- 
tice displayed in fig. |2j which are identical everywhere except 
on one face. Because the ground state (Hli is written in terms 
of chromatic polynomials, the identity (|9]l results in a local re- 
lation between amplitudes in the ground-state wavefunction: 



(10) 



The extra factor of y/cj) in jiojl relative to (jojl results from the 
weight per unit length in (|4k!^ 

I I ~ J L 



FIG. 2: Three configurations mixed by the Jones-Wenzl projector 



We can exploit the relation (lOi to find a projector Hj 
that annihilates the ground state 



loops. It is of the form Hj — 



W) and couples the winding 
Jout, and acts non-diagonally 



on the four "inside" links around the square in fig. |2] and 
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diagonally on the remaining four "outside" links. The non- 
diagonal part of each Hj acting on the inner links is chosen so 
as to project onto the linear combination (t>^^'^\t) — 4>\E) + |/). 
A Hermitian operator doing this is 



(11) 



acting on \t), \E) and |/) respectively; it annihilates all other 
configurations on the inner links. The diagonal part Jout an- 
nihilates any configuration on the outer links other than those 
illustrated in fig.|3] i.e. labeling the outside links ji, j2, js and 
j4 from left to right and letting Qj — {1 — Pj ) gives 

Jout — QjlPj2Qj3J^ji + Qj2J^jiQ j^Pji 



FIG. 3: Configurations on the outside links not annihilated by Jout 

Each Hj therefore mixes each of four sets of three config- 
urations among themselves; one of of these sets is illustrated 
in fig.[2j while the other three sets are given by changing the 
configurations on the outer links to one of the others in fig. 
[3] Because of the presence of Jout, Hj does not mix any net 
configurations with any non-net configurations. Using (10 1 



then ensures that Hj\^>) = for j any set of eight links of the 
form illustrated in fig. [2] i.e. the four spins on a face and the 
four touching two adjacent vertices. There are thus four such 
terms for each face on the lattice. Since j^*) can be rewritten 
in terms of nets on the dual lattice as well, the analogous pro- 
jectors H-j also anirihilating j^*) can be defined by repeating 

the above arguments acting on the |0);, basis on the dual 
lattice. The full Hamiltonian is then 



(13) 



where e is a coupling strength included for convenience. 

Numerical results — We have confirmed by exact diago- 
nalization that the Hamiltonian ( pj) has exactly four ground 
states when e > 0, as required for doubled Fibonacci topo- 
logical order. More importantly, we have also checked that 
H is gapped as for e > 0. Our data is shown in figs. |4]and 
|5] In fig. |4] we show the gap A to the first excited state as a 
function of e for different system sizes. To illustrate how the 
gap survives in the thermodynamic limit, we show in fig.|5]its 
finite-size scaling. This suggests that the gap is indeed non- 
zero, with an estimate of A w 0.025 in the thermodynamic 
limit for e = 0.8. We checked as well the gap including the 
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FIG. 4: The gap to the first excited state as a function e for different 
system sizes. 
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FIG. 5: The gap to the first excited state as a function of the inverse 
system size at e = (circles, no Jones-Wenzl terms) and at e = 0.8 
(squares). Lines guide the eye. The system is gapless at e = and 
gapped at e = 0.8. 



Hj but omitting the dual projectors Jfj, and obtained about 
half this value. 

We diagonalized H by the Lanczos method, using a paral- 
lel and highly scalable exact diagonalization code [21]. The 
spins live on the links of a Li x L2 square lattice with peri- 
odic boundary conditions, so that the clusters contain 2L1L2 
sites. Clusters with periodic boundary conditions of size 3x3, 
VTO X yiO, 4x3, %/l3 X Vl3, 5 X 3, and 4 X 4 were di- 
agonalized. The 3x3, VW x v^, ^/T3 x ^13, and 4 x 4 
clusters have the symmetry of the infinite lattice. We exploit 
translation symmetry to reduce the Hilbert-space size, so the 
largest Hilbert space size is about 2.7 x 10^ for the 4x4 lat- 
tice. It is not very large, but we have complicated multi-spin 
interactions and there are many matrix elements per row; e.g. 
roughly 16000 for the 4x4 lattice. The lowest eigenvalues 
for e > are in the zero-momentum sector 
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Conclusion- We have presented strong evidence for a 
magnetic model with a gap and doubled Fibonacci topolog- 
ical order Our results illustrate how to simplify the string-net 
model while still remaining in the same topological phase. We 
believe this is the "minimal" gapped Hamiltonian with an ex- 
act ground state in this phase, thus completing the program of 
Ref 8 in the Fibonacci case. By expressing the ground state in 
terms of a topological quantity, the chromatic polynomial, our 
construction also illuminates precisely how topology appears 
in this ground state. Moreover, since H is the sum of non- 
commuting projectors, the model is not exactly solvable, so 
that the excitations are dynamical, and not just fixed defects. 

It would be very interesting to understand how this Hamil- 
tonian could be simplified further (e.g. by omitting Jout in 
Hj) while still preserving the same topological order It would 
also be illuminating to apply a similar analysis to other mod- 
els with non-abelian topological order. A good starting point 
may be the models of e.g. Refs.|22l[23] related to the one con- 
sidered here. Another interesting topic for further study is the 
presumed quantum critical point at e = describing a tran- 
sition out of the topological phase; an analogous point also 
occurs in the correspondingly deformed toric code. 

The work of RF. is supported by the US National Science 
Foundation under the grant DMR/MPS 1006549, while the 
work of S.I. and M.T. is supported by the Swiss HP2C initia- 
tive. M.T. acknowledges hospitality of the Aspen Center for 
Physics, supported by NSF grant # 1066293. The calculations 
were performed with the MAQUIS -ED code, developed with 
support of the HP2C, on the Brutus cluster at ETH Zurich and 
the Cray XE6 "Monte Rosa" at CSCS. 



[1] A .Y. Kitaev, Annals of Phys. 303, 2 (2003) |arXiv:quant-] 
[~I)h/9707021| 
[2] C. Nayak, S. H. Simon, A. Stem, M. Freedman and S. Das 
Sai-ma, Rev. Mod. Phys. 80, 1083 (2008) |arXiv:0707.1889 | 



[3] M.A. Levin and X.-G. Wen, Phys. Rev. B 71, 045110 (2005) 

|arXiv:cond-mat/0404617| 
[4] I. Klich, Ann. Phys. 325, 2120 (2010) |arXiv:0912.0945| 
[5] S. Bravyi, M. B. Hastings, and S. Michalakis, J. Math. Phys. 51 

093512 (2010); S. Bravyi and M. B. Hastings, Comm. Math. 

Phys. 307 , 609 (2011) |arXiv: 1001.43631 
[6] P Fendley, Ann. Phys. 323, 3113 (2008) | ar Xiv:0804.0625 J 
[7] P Fendley and J. L. Jacobsen, J. Phys. A 41, 215001 (2008) 

|arXiv:0803.2618|. 
[8] M. Freedman, Comm. Math. Phys. 234, 129 (2003) 

I arXiv:quant-ph/01 10060 1 
[9] M. Freedman, C. Nayak and K. Shtengel, Phys. Rev. Lett. 

94, 066401 (2005) |arXiv:cond-mat/0312273|; ibid9A, 147205 

(2005) |arXiv:cond-mat/0408257|; Phys. Rev. B 78, 174411 

(2008) |arXiv:cond-mat/0508508| 

[10] S. Trebst, P Werner, M. Troyer, K. Shtengel and C. Nayak, 

Phys. Rev. Lett. 98, 070602 (2007) | arXiv:cond-mat/0609048 1 
[11] M. Troyer, S. Trebst, K. Shtengel and C. Nayak, Phys. Rev. 

Lett. 101, 230401 (2008) |arXiv:0805.2177|. 
[12] D. Rokhsar and S. Kivelson, Phys. Rev. Lett. 61, 2376 (1988) 
[13] R. Moessner and S.L. Sondhi, Phys. Rev. Lett. 86, 1881 (2001) 

I arXiv:cond-mat/0007378 1 
[14] M. B. Hastings, Phys. Rev. B 69, 104431 (2004) [arXiv: corid^ 

mat/0305505 1; Europhys. Lett. 70, 824 (2005) [arXiv: cond^ 
[_ mat/04 11 094 1 

[15] P Fendley and E. Fradkin, Phys. Rev. B 72, 024412 (2005) 

I arXiv:cond-mat/050207 1 1 
[16] L. Fidkowski, M. Freedman, C. Nayak, K. Walker and Z. 

Wang, Comm. Math. Phys. 287, 805 (2009) lai'Xivxond^ 

mat/0610583 1 

[17] FJ. Burnell and S.H. Simon, Annals Phys. 325, 2550 (2010) 

I arXiv: 1004.55861. 
[18] P Fendley a nd V. Krush kal, Geometry & Topology 13, 709 

(2009) I arXiv: 071 1.00161 

[19] V.FR. Jones, Invent. Math. 72, 1 (1983); H. Wenzl, 
C. R. Math. Rep. Acad. Sci. Canada 9, 5 (1987) 

[20] W.T Tutte, Graph theory as I have known it (Oxford University 
Press, New York, 1998) 

[21] S. V. Isakov and M. Troyer, to be published. 

[22] A. Velenich, C. Chamon, X.-G. Wen, J. Phys. A 43, 172002 

(2010) |arXiv:0911.1132| 
[23] B. Paredes, arXiv: 1209.6040, 



